QUANTIZATION OF SINGULAR SYSTEMS AND 
INCOMPLETE MOTIONS 



N.P. LANDSMAN 

Korieweg-de Vries Institute for Mathematics, University of Amsterdam, Plantage 
Mmdergracht 24, 1018 TV AMSTERDAM, THE NETHERLANDS 
E-mail: npl@wins.uva.nl 

The need for a mathematically rigorous quantization procedure of singular spaces 
and incomplete motions is pointed out in connection with quantum cosmology. We 
put our previous suggestion for such a procedure, based on the theory of induced 
representations of C*-algebras, in the light of L. Schwartz' theory of Hilbert sub- 
spaces. This turns out to account for the freedom in the induction procedure, at 
the same time providing a basis for generalized eigenfunction expansions pertinent 
to the needs of quantum cosmology. Reinforcing our previous proposal for the 
wave-function of the Universe, we are now able to add a concrete prescription for 
its calculation. 



1 Introduction 

1.1 Incomplete dynamical systems 

Minisuperspace cosmologies with singularities are examples of dynamical sys- 
tems whose motion may be incomplete (that is, in a given parametrization the 
flow may not be defined for all i£K). In general, classical incomplete motion 
may occur even on smooth phase spaces; for example, a particle may escape to 
infinity in a finite time. In cosmological examples, though, the incompleteness 
is generically caused by the presence of singularities in the space of physical 
degrees of freedom of the model in question. The notion of singularity used 
here is rather wide; boundary points to a manifold are included as possible 
singular points. Indeed, the incompleteness of the motion of a cricket ball off 
a batsman hitting a boundary is precisely caused by the boundary of the field. 
On the other hand, one may think of the incompleteness of the motion of an 
observer falling into a black hole. The question whether the space is question 
may be extended so as to render the given motion complete is interesting, but 
not quite relevant in this context. 

In the general theory of dynamical systems, incompleteness is not much 
studied, for the simple reason that under mild conditions one may reparametrize 
the flow, so as to make the motion complete. The question whether a singu- 
larity is approached in finite time therefore appears to be ill-defined. However, 
the physics of the situation usually singles out a preferred choice of the time 
co-ordinate, an issue that has been much discussed in classical and quantum 
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gravity; see IshamEJ and references therein. 

In our opinion, in a canonical context the choice of a time parameter is 
determined by the choice of the physical Hamiltonian, which simply defines 
time through Hamilton's equations. Hence the above question in any case 
makes sense in the setting of Hamiltonian dynamical systems. Further to the 
well-studied classes of integrable and chaotic systems, one should therefore 
study the class of incomplete Hamiltonian dynamical systems. While this is 
fascinating already classically, the issue of genuine physical importance lies 
in the quantum theory of such systems. It is well known that, in view of 
the unitarity of the time-evolution, a quantum system can be neither chaotic 
(in the classical sense) nor incomplete. Hence one faces the difficult task of 
identifying properties of the quantum Hamiltonian which indicate that the 
classical limit of the theory is integrable, chaotic, or incomplete. Signatures 
for chaos are found in the statistical distribution of the eigenvalues, whereas 
classical incompleteness is often related to the fact that the Hamiltonian fails to 
be essentially self-adjoint on its natural domain of smooth compactly supported 
wave- functions 0. It is intuitively evident why this is so: such wave- functions 
have support away from the singularities, so that the natural domain contains 
no information about the boundary conditionsB. 

1.2 Singular reduction 

Our second source of inspiration lies in the origin of the singularities of the 
physical phase space of general relativity. Namely, since the space of unphysical 
degrees of freedom (in whatever formalism one uses) is regular, the singular- 
ities come from the constraints of Einstein gravity. The singularities of the 
reduced phase space are well understood: the physical phase space is strati- 
fied by symplectic manifolds, each of which is stable under the Hamiltonian 
flow generated by any function that Poisson-commutes with all constraints □. 
This important infinite-dimensional fact is entirely analogous to the case of 
singular Marsdcn-Wcinstcin quotients by proper actions of finite-dimensional 
Lie groups on finite-dimensional symplectic manifolds 

When the symplectic space one starts from is a cotangent bundle S = T*Q, 
and the group action on S (with equivariant momentum map J) is pulled 
back from a G-action on Q, the Marsden-Weinstein quotient at zero is simply 
J _1 (0)/G = T*(Q/G), where away from the singular points the cotangent 
bundle is defined as usual, whereas at the singular points it is defined by 
the left-hand side. One may quantize the reduced space by saying that the 
pertinent Hilbert space is L 2 (Q/G), and that the Hamiltonian is minus the 
Laplacian, defined on the domain C%°(Q/G)B. As already mentioned in the 
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previous subsection, this operator is not essentially self- adjoint; one needs to 
choose boundary conditions, so there is no unique quantum theory. 

Instead, one needs a method of constrained quantization which parallels 
the classical procedure of symplectic reduction, and which still applies when 
the usual conditions guaranteeing that the redjured space be smooth □ do not 
apply. The method proposed by the authorMtJ satisfies these requirements, 
and simultaneously solves the problem described in the previous paragraph. 

The essence of this method is best explained bjz-comparing it with Dirac's 
well-known technique of constrained quantizationl!3. We restrict ourselves to 
first-class constraints. Recall that symplectic reduction is a two-step proce- 
dure: firstly, the constraints ipi — are imposed, and secondly, roughly speak- 
ing, gauge-equivalent points on the constraint hypersurface are identified. In 
quantum theory only one of these steps has to be taken. Dirac's approach 
to constrained quantization selects the first step, in imposing the quantized 
constraints (pi as state conditions (p^ — on the Hilbert space H of the 
unconstrained system. 

These equations rarely have solutions in 7i, and are, accordingly, usually 
solved in some enlargement V of TL 13. Since the inner product of TL is not 
defined on V, this leads to certain problems. Moreover, it is not a priori clear 
which enlargement to choose. For example, when ipi is a second-order differ- 
ential operator, a case which is of prime importance for quantum cosmology, 
there will be two linearly independent eigenfunctions \&q, i = 1,2, with eigen- 
value zero. In cosmological models, where p>i is the Wheeler-DeWitt operator 
(that is, the quantized super-Hamiltonian), the wave- function of the Universe 
is supposed to be a certain linear combination of and The formalism 
doesn't tell which combination to choose. 

Partly in response to these problems, and partly on the basis of purely 
mathematical considerations, the author's method of constrained quantization 
singles out the second step of the classical symplectic reduction procedure as 
the one to be quantized. We will first review this method in its original for- 
mulation, and then introduce a refinement. Upon the latter, our method turns 
out to be surprisingly closely related to Dirac's method, applied to the enlarge- 
ment V, without sharing the difficulties just mentioned. A related approach is 
discussed by Marolf fcHl. 

2 Quantized symplectic reduction 

Applied to the special case of Marsden-Weinstein reduction of cotangent bun- 
dles, as above, our method amounts to choosing a dense subspace £ of L 2 (Q) 
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on which the quadratic form 

(*,*)o= / dx(%U(x)$)v (Q) (1) 
Jg 

is well defined. Here U is a unitary representation of G on L 2 (Q) which 
quantizes the G-action on Q, and dx is a Haar measure on G. Unless G is 
compact, this form tends to be unbounded and even non-closable on £, but 
one can nonetheless form the quotient £/Af of £ by the null space AT of ( , )o- 
When the latter is positive semi-definite, £jM is a pre-Hilbert space in the 
sesquilinear form inherited from ( , ) . Its completion 7i° is the physical state 
space of the constrained system, and as such serves as the quantization of 
classical reduced space J _1 (0)/G = T*(Q/G). This is true whether or not the 
Marsden-Weinstein quotient is regular. 

The classical physical observables are functions F on T*Q which Poisson- 
commute with the G-action. The quantization F of such an observable on 
the unphysical Hilbert space L 2 (Q) should commute with U(G), and leave £ 
stable. If these two conditions hold, F leaves Af stable, and therefore quotients 
to an operator F° r pn £/Af. If F is bounded, F° is bounded under suitable 
assumptions on Gu, and may be extended to Tt° by continuity. When F is 
unbounded, containing £ in its domain, this procedure yields an unbounded 
operator on the dense subspace £/Af of H°. The point is now that in typical 
examples F° tends to be essentially self-adjoint on the projection £ of £ to 
H°. This is possible, because £° is generally larger than C%°(Q/G). ._. 

To see this in a simple example, classically due to Gotay and Bos E3, 
consider the standard action of G = SO (2) on Q = M 2 . The momentum map 
for the pull-back action on S — T*M. 2 is J(p, q) = q A p, so that the quotient 
J- 1 {0)/SO(2) may be identified with x E (where M£ := E+ U 0), which 
may be thought of as the cotangent bundle T*Rf • Hence the singularity in 
the reduced space takes the form of a boundary. 

Quantum reduction is done with H = L 2 (R 2 ) and £ = G^°(M 2 ). The 
unconstrained quantum Hamiltonian H phyB = —A + V(r) on L 2 (M. 2 ) is SO(2)- 
invariant and essentially self-adjoint on G^°(1R 2 ). Using the unitary transfor- 
mation U : L 2 (R + ,rdr) — > L 2 (R + ,dr) defined by U~$>(r) := ^Jr^(r), we have 

UHL»U* = --^-^+V(r). (2) 

While the analysis of this expression is quite straightforward for any reasonable 
potential V, the free case V = already suffices to illustrate the main point. 

Defined on G C °°(R+), the operator (|) is in the limit circle case el at and 
in the limit point caseel at oo. Hence it has deficiency indices (1, 1), so that it is 



4 



not essentially self-adjoint. However, defined on £° — C£°(IR 2 ) , which consists 
of functions of the type \t(r) = y/rf(r 2 ) with / 6 C^°(M.q), the operatorJxi 
question is essentially self-adjoint, as follows from Thm. 3 in Nussbaum 
The closure of the latter operator is an extension of the closure of the former, 
to whose domain one adds functions of the indicated type in order to achieve 
essential self-adjointness. The boundary condition ^(0) = corresponds to a 
hard wall potential at the origin. 

Other examples of this phenomenon are given by WrentLj, who looks at 
the quantization of Stieffel chambers (i.e., quotients of a maximal torus of 
a compact Lie group by its Weyl group), finding that our procedure assigns 
Neumann boundary conditions to the Laplacian. 

3 Minisuperspace quantum cosmology 

Let us now apply our quantization method to the simplest cosmological models 
in minisuperspace. Using the canonical formalism (see our earlier work El for 
a covariant approach), we assume that the Hamiltonian constraint H = 0, 
in which H is a scalar, is the only constraint on the classical minisuper phase 
space. If the Hilbert space TL results from the quantization of the unconstrained 
phase space, and the Wheeler-DeWitt operator H is the quantization of the 
Hamiltonian constraint (realized as an unbounded self-adjoint operator on TL), 
our formalism stipulates that firstly one should find a dense subspace £ C TL 
on which the quadratic form 

(tf,$)o = / dt(^>,e UH <$>) (3) 
Jr 

(cf. (|l|)) is well-defined. Secondly, one should determine the reduced Hilbert 
space TL°, defined as in the preceding section. It would seem, therefore, that 
one needs to explicitly compute the null space J\f of (, )o- Fortunately, in 
practice there is often no need to do so. One may, instead, start from an 
Ansatz Hilbert space TL°, and find a 'quantum reduction map' V : £ — > TL° 
satisfying 

(V%V*) iiP = (%$) . (4) 

It is easily seen that V quotients and extends to a unitary isomorphism between 
TL° and TC°, through which one may transfer the physical Hamiltonian and 
other observables from the unknown space TL° to TL , where everything is 
explicitp. The operators A on TL that commute with H and leave £ stable, are 
represented by physical observables A on TiP by the prescription 

A°VV = Ui°*. (5) 
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In the case at hand, these two steps. may be taken at one stroke. Using 
the theory of eigenfunction expansionstrlLj, one finds a Hilbert space Ti- and a 
Hilbert-Schmidt injection Ti <—* Ti- , such that (almost) all relevant generalized 
eigenfunctions of H (that is, those contributing to the spectral resolution of 
) lie in Ti-. Hence Ti is a Hilbert subspace of Ti- in the sense of Schwartz 
(i.e., the embedding Ti Ti- is continuous). If Ti+ is the continuous dual 
of Ti-, one obtains an anti- linear map W — > "J from Ti + to Ti, defined by 

(*,$) = 0M) W (6) 

for all $ £ Ti. Here the restriction of W £ Ti+ to Ti C Ti- (with its own Hilbert 
space topology) is continuous because of the continuity of the embedding Ti «-> 
Ti- ; the Riesz-Fischer theorem then implies the existence of an element ^ £ Ti 
for which (g) holds. The image Ti+ of Ti + in Ti under this map is always dense, 
and when Ti is dense in Ti-, the map \& — > $ is injective, so that 7i + may 
be seen as a dense Hilbert subspace of Ti. One then has a Gel'fand triplet (or 
rigged Hilbert space) Ti+ C Ti C Ti-. However, it is by no means necessary 
that Ti be dense in Ti- , so that the formalism of Schwartz is a generalization 
of that of Gel'fand. 

In case that the spectrum cr(H) does not contain 0, the physical Hilbert 
space Ti° determined by (||) is empty. When H has in its discrete spectrum 
(Jd{Ti), so that the Wheeler-DeWitt equation H^> = has a solution ^ £ 
Ti, there is no need for our formalism. This rarely occurs in minisuperspace 
models; instead, we assume that lies in the absolutely continuous spectrum 
fac(W). We may take Ti° as the corresponding multiplicity space, realized in 
some arbitrary fashion, with 

V*(m) = <*,VV(0)>, (?) 

defined for \P £ Ti+, so that f e£. Here y M (A), where A £ u(E) and fi labels 
a basis of Ti°, is a generalized eigenfunction of H in Ti- with generalized 
eigenvalue A. Note that ^ is well defined, for the <p M (A) lie in the closure of Ti 
in Ti-ft One easily verifies (0) from ([|) and eq. (1.17) in §V.l of Berezanskii's 
bookB. 

An operator A on Ti that leaves Ti+ stable has a dual A* : Ti- — > W_ . If, 
in addition, A commutes with H , the dual A* cannot change the generalized 
eigenvalue A, so that, in particular, one has A*(^ M (0) = A* for certain 

coefficients A*^ v . Combining this with ([?]) and (|B|), we see that an observable 

A phyB on Ti is realized on the physical state space Ti° by 

i° hys /(M) = 3^/M- (8) 
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4 The wave-function of the Universe 

Physical processes in the quantum Universe may be computed if one specifies 
the physical Hamiltonian; as pointed out before, this should be an operator 
H phyB on the unconstrained Hilbert space Ti which commutes with the Wheeler- 
DeWitt operator H, and leaves Ti+ stable. This operator induces an operator 
H® hya on the physical Hilbert space Ti , given by (J|). One may anticipate that 
the spectrum of H° hys is purely continuous, since otherwise the Universe would 
settle in its ground state, and nothing would happen. 

For example, a classical homogeneous and isotropic universe filled with 
noninteracting dust may be described by the phase space T*R 2 , with co- 
ordinates (a, 4>) standing for the logarithm of the radius and the dust field, re- 
spectively, with conjugate momenta (pa,p<f>)- The classical super-Hamiltonian 
is 

H K = ±{pi + K e 4a -pl), (9) 

where k = or ±1. A reasonable choice for the classical physical Hamiltonian 
is 

H phyB = pf, (10) 

the idea is that the dust field serves as a physical time parameter^, but since 
4> itself does not Poisson-commute with the constraint H K , this idea must be 
implemented through its conjugate p^, which does. 

The unconstrained Hilbert space is simply Ti = L 2 (R 2 ), with wave- functions 
^(a,4>). The Wheeler-DeWitt equation 

i (-{d/daf + Kexp(4a) + (d/dc/)) 2 ) * = (11) 

can be solved explicitly for all values of tM; we now add the knowledgeJJiat 
the relevant generalized eigenfunctions must be polynomially bounded E3 in 
order to contribute to the spectral resolution of H. For k = this leads to 
Hi = L 2 (R) <g> C 2 , whereas for k = ±1 one has U° ±1 = L 2 (R). The generalized 
eigenfunctions at are, omitting the argument (0), 

pgW) =exp[ifc(a± «£)]; (12) 
ifl^ia, 4>) = ^- 1 e^Vsinh(^|fc|/2)X l | fe | /2 (ie 2Q ); (13) 
¥>r -1 («. 4>) = ie^Vcosech(7r|fe|/2)(J i | fe | /2 + J_ J | fc |/ 2 )(ie 2 «). (14) 

With (|To| ) as the physical Hamiltonian, it is reasonable to ask the wave- 
function of the Universe to peak around k — 0; this approximately yields 
<f K=0 ~ 1, * K=1 ~ K , and * K=0 ~ J . This agrees with Hartle-Hawking 
boundary conditions for k = —1 only, and always disagrees with Vilenkin 
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boundary conditions (as imposed by ZhukES). However, a direct comparison 
is difficult, because these authors use a mathematical framework unrelated to 
the theory of self-adjoint operators on a Hilbert space. 
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